arXiv:hep-lat/0503014v2 18 Jul 2005 


Finite volume effects 
for meson masses and decay constants 


Gilberto Colangelo, Stephan Diirr and Christoph Haefeli 

Institut fiir Theoretische Physik, Universitat Bern, Sidlerstr. 5, 3012 Bern, Switzerland 

14 March 2005 


Abstract 

We present a detailed numerical study of finite volume effects for masses and decay 
constants of the octet of pseudoscalar mesons. For this analysis we use chiral perturbation 
theory and asymptotic formulae a la Liischer and propose an extension of the latter beyond 
the leading exponential term. We argue that such a formula, which is exact at the one-loop 
level, gives the numerically dominant part at two loops and beyond. Finally, we discuss 
the possibility to determine low energy constants from the finite volume dependence of 
masses and decay constants. 
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1 Introduction 


In lattice QCD the determination of the mass and decay constant of the lowest-lying state 
with a given set of qnantnm nnmbers is entering the high-precision era. Even in the fully 
unquenched case (i.e. with sea- and valence-quarks being degenerate) the pion mass can be 
measured, for hxed bare parameters, with an accuracy at the percent level, and future progress 
towards the permille level is anticipated. However, to make contact with the real world, three 
extrapolations are needed. These are {i) the continuum extrapolation, (ii) the inhnite volume 
extrapolation, and {Hi) the chiral extrapolation. In each case there is considerable help and an 
analytical guideline from an effective held theory framework. 

This article is concerned with the extrapolation to inhnite volume, where the situation is 
particularly favorable (for a compact review of the recent literature see P). As shown by 
Gasser and Leutwyler, chiral symmetry imposes strong constraints on the dynamics at low 
energy in QCD, even if the system is enclosed in a hnite box. Accordingly, Chiral Perturbation 
Theory (ChPT) may be adapted to the hnite volume case |2l El El- In this framework hnite 
volume ehects can be taken systematically into account, in a perturbative loop expansion. If 
the spatial volume is large enough internal degrees of freedom of the particle of interest play 
no role as far as hnite volume ehects are concerned and these are exclusively due to pion loops. 
This means that they hrst appear at next-to-leading order (NLO), i.e. at O(p^) in the chiral 
counting. Another consequence is that they are exponentially small in the pion mass for any 
particle that couples to the pion held: the ehect behaves like for pions, kaons, etas and 

nuclei. To date, a number of hnite volume calculations have been performed at one-loop order. 
In the nucleon sector the corrections to the mass iniini, the magnetic moment [Hj and the baryon 
axial charge jZj have been worked out. In the meson sector the original calculation for Mjr(L) 
and T)r(T) |2] has been extended to the quenched case jH] . The same quantities have also been 
analyzed in a quark-meson model [H]. More recently, the hnite volume shifts for Fk{L) and 
Bk{L) have been given in Ref. pH] and the extension to heavy-meson chiral perturbation theory 
has been described m- No full two-loop calculation of hnite volume ehects has appeared yet. 

For some observables the Liischer formula represents a convenient and powerful alternative 
jl2j . It allows to estimate subleading (in the chiral counting) hnite volume ehects for the mass 
of a particle with less ehort, while sticking to the leading order in an expansion in powers of 
g-Mvri fppg formula gives the hnite volume shift Mp{L) — Mp of a particle P in terms of the 
infinite volume ttP forward scattering amplitude in the unphysical (Euclidean) region. For 
this amplitude the ChPT expression at a certain loop order is used. The approach via the 
Liischer formula is economical, since only a chiral calculation in inhnite volume is needed, and 
the loop in hnite volume comes for free. Applications of the Luscher formula to the mass of the 
nucleon 131131 and the pion ^31 have been worked out. In the latter publication several orders 
in the ChPT input for the tttt scattering amplitude were compared and it was found that in a 
certain range of and L subleading ehects (in the chiral counting) can be large with respect 
to the leading contributions. An extension of the Luscher formula has been constructed for 
pseudoscalar decay constants ng. 

This article presents a resummed version of the Luscher formulae for masses and decay 
constants (this has been briehy discussed by one of us in [T]), where the terms neglected in the 
large-L expansion are ^ rather than We proceed with the asymptotic 

expression for Mt,{L)— to 3-loop order, and for Ft,{L) — Mk{L) — M^, Fk{L) — Fk to 
2-loop order, by using the available knowledge in the literature on the scattering amplitude or 
axial-vector matrix element. In all cases, the result may be given in a compact formula that 
does not involve any numerical integration. 
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2 Finite volume effects 


Lattice calculations are necessarily done in a finite 4D volume which acts as an IR-cutoff. 
Typically a. L^xT geometry with periodic boundary conditions in all directions is chosen with 
T ^ L and both large compared to the inverse temperature of the QCD phase transition or 
crossover. In order to determine the mass of a particle one considers a correlator of two properly 
chosen interpolating helds 

/ CXD 

d^x {(j){x)(j){0)) (1) 

n=0 

and tries to determine the energy levels (typically the lowest) at zero spatial momentum, 
Mn{L) = En{p = 0). In the following, we will be concerned with the shift M{L) — M, where 
M = M{L = oo), for the groundstate (n = 0) due to the finite 3D volume L^. Decay constants 
or, more generally, matrix elements suffer from analogous shifts at finite spatial volume. For 
not-too-small box length L these shifts can be calculated analytically, thus offering a means 
to correct lattice data for this systematic effect. We now give a brief outline of the two main 
frameworks for such a calculation. A comment on cut-off effects and their interplay with hnite 
volume effects is given in appendix iBl 


2.1 ChPT in finite volume 


In QCD with light flavors the physics in the infrared region is controlled by chiral symmetry. 
As shown by Gasser and Leutwyler this still holds true if the system is enclosed in a hnite box 
xT, provided both L and T are large enough that chiral symmetry is not restored j3]. They 
have shown that in an isotropic box with periodic boundary conditions for the meson helds the 
hnite-volume dependence comes in exclusively through the propagators. The latter becomes 
periodic in all spatial directions, and in the limit T ^ oo can be written as follows 

G(x°,x) = ^Go(a:°,x-h nL) (2) 


which is equivalent to replacing the integration over the spatial part of the momenta by a sum 
over multiples of 27i{L. In other words, with periodic boundary conditions the Lagrangian 
remains the same as in inhnite volume. 

The expansion parameters in ChPT are 


dvriQ ’ dvriQ 

and the theory can be meaningfully applied only if both are small. In a hnite volume spatial 
momenta are discretized: p = 27rn/L with n a vector of integers. Therefore one can have 
“small” nonzero momenta and apply ChPT only if the condition 
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~ 1 fm 


(4) 


is satished. A priori there is no way to say how much L has to be in excess of 1 fm. As a 
guideline we observe that the lowest non-trivial momentum in a 1 fm box is 1.2 GeV, which 
is certainly beyond the realm of ChPT. Note, hnally, that unlike iQT the combination 
is not constrained. Both -C 1 and 3> 1 are acceptable n El ii cni, but they imply 
diherent ways to organize the chiral series. 


MttT 3> 1 ^ “p-expansion” 


e-expansion 


2 


(5) 

( 6 ) 





n 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

m{n) 

6 

12 

8 

6 

24 

24 

0 

12 

30 

24 

24 

8 

24 

48 

0 

6 

48 

36 

24 

24 


Table 1: The multiplicities m{n) in (I12j) for n < 20. 


Here we shall restrict ourselves to the former case, where the chiral counting is 

Ml m 0{p‘^), 1/L P O (p) . 


(7) 


With this setup Gasser and Leutwyler calculated the mass and decay constant shift in a theory 
with — ^ degenerate pseudo-Goldstone bosons, and obtained |2j 


M^(L) 

F^{L) 


F. 


Nf 

1“ + 0{i^) 


Here we have introduced the abbreviations (note the in the denominator for all P) 

^ Ml 
- (4vrF.)2 

\p = MpL 

for P = 'K (and P = K^r] will be used below) as well as the modified ^ shape function 

^ 4 m(n) ^ , 

9i{x) = 2^ Kii^nx) 


n=l 


nx 


( 8 ) 

( 9 ) 


( 10 ) 

( 11 ) 


( 12 ) 


where Ki is a Bessel function of the second kind and the multiplicities m{n) have been given 
in ^3], but for convenience we reproduce them in Tab.^ Given the asymptotic expansion 
Ki{z) ~ tt/ {2z)e it is clear that in the p-regime eqns. fjH t l^ fT^ represent quickly converging 
expressions. Several observables have been worked out at one-loop order iniininoiiniMEi 
[TH], but to date no two-loop result obtained in this setup has appeared. 


2.2 Liischer formula 

An entirely independent approach has been devised by Liischer who has proven an elegant 
relation between the mass shift of the particle P in a hnite volume L^xcx) and the Pn scattering 
amplitude in inhnite volume US 

Mp(L) - Mp = --4^ r dyTp(iy)e-'/^'^ + 0(e-^‘‘) . (13) 

lOTT^Ap J-oo 

Here Fp{i') denotes the inhnite volume forward {t = 0) scattering amplitude of P and tt in 
Minkowski space. The integration runs along the imaginary axis, i.e. Fp{v) is evaluated for 

u = iy (14) 

with real p, thus staying far away from the cuts. Only the real part of Fpiiy) contributes to 
the integral, since the imaginary part is odd in y. An additional piece in the original formula, 

^Our gi relates to gi of [2] via giiM) = (47r/M,r)^-(7i(M,r,/3 = oo, L) and is a dimensionless function. 
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referring to the 3-particle vertex, is omitted here, since we assume P to be a pseudo-Goldstone 
boson. The Liischer formula keeps only the leading term in an expansion in (fractional) 
powers of The generic bound M > 2Mn- can be specihed to M = in a theory 

with pseudo-Goldstone bosons only. 

Recently, an analogous “Liischer-type” formula has been derived for the hnite volume shift 
of the axial-vector decay constant Fp. It reads na 

Fp(L) - Fp = +— Ir r ‘‘y e-v'^^y + 0(e-"'-) (15) 

oTT^Xp J-oo 

where A/p(z/) is derived from the matrix element (yrTrl^d^lP) via a subtraction prescription we 
will specify below. Like in the mass formula 03 the hnite volume shift of Pp is expressed in 
terms of an inhnite-volume amplitude, evaluated in the unphysical (Euclidean) region, thus far 
away from the cuts. Again, only the leading term in an expansion in (fractional) powers of 
is kept. Note the reverse overall sign, compared to (d, and the fact that the net physical 
effect is opposite to this prefactor; in other words Mp(L) > Mp and Fp{L) <Fp. 

To predict the shifts Mp{L)—Mp and Fp{L)—Fp in a lattice calculation with a known box 
length L, and thus to correct the data for this systematic effect, the formulae dd must be 
fed with an explicit representation of the amplitudes Fp{v) and N'p{i>), respectively. This is the 
place where GhPT naturally enters, even if one opts for the Liischer approach. Using existing 
knowledge about the relevant amplitude at n-loop order, one gets the leading piece, in the 
expansion, of the hnite-volume shift of Mp{L), Fp{L) to n-|-1-loop order. For instance, using 
the tree-level expressions Ftp{v) = —MI/F^ and = —2Mtp/F^^ in 2-flavor GhPT yields 

3 /\/f2 

M,(L)-Mp = + ( 16 ) 

F,(L)-F, = 

in agreement |2j with the 1-loop chiral expressions (jH)) and ©. Because of the “elevator”-effect 
in the loop expansion and since the associated chiral calculation is in inhnite volume, it is much 
easier to push to higher chiral orders in the Liischer-type setup than with a straightforward 
GhPT-in-£nite-volume calculation |Tn]. This offers a genuine opportunity to compare several 
chiral orders and thus to assess the chiral convergence behavior. Indeed, in m the hnite 
volume shift in the pion mass was evaluated, using GhPT at LO/NLO/NNLO for Ftp{v) to get 
the asymptotic piece of the full chiral expression at 1/2/3-loop order, and it was found that 
for some (M^, L)-combinations the chiral series converges well, if at least the NLO input is 
included. Still, one might worry whether the non-asymptotic pieces of order 0(e“^^’") would 
prove numerically relevant • In ^Ij a hrst attempt was made to discriminate those regions in 
the (Mtt, L)-plane where higher orders in the chiral expansion dominate against those regions 
where terms omitted in the Liischer approach are more important. Below, we shall present 
a resummed version of the Liischer-type formulae (HI ESI), where the pieces cx e and 
g^g-v^A,r included and the terms estimated. On this basis a more precise 

assessment of the relevance of higher loop corrections versus higher powers of can be made. 

Note hnally that all Liischer-type formulae build on the unitarity of the theory and thus 
hold for the full (unquenched) theory. In the (partially) quenched case it seems indispensable 
to start in the framework of subsect. 2.1, but even then the arguments for using the inhnite 
volume Lagrangian reside on less solid grounds - see Ref. [H] for a lucid discussion. 
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3 The Liischer formula resummed 


In Liischer’s derivation of the asymptotic formula for the hnite volume correction to particle 
masses the first step is a proof that the leading exponential term is given by the sum of all 
diagrams in which only one propagator is taken in hnite volume. This class of diagrams yields 

Mp{L) -Mp = Y. J G'o(g)r(p, q, -p, -q) + ... (18) 

where Go{q) = 1/{M^ + q^) is the full propagator and T{j)i,p 2 ,Pz,Pi) the four-point vertex 
function in inhnite volume. Liischer then concentrates on the leading exponential contributions 
(those with |n| = 1), and shows that, if one disregards terms which are exponentially suppressed 
with respect to exp(—A t^), three of the four integrations in (fT^ can be performed explicitly and 
the result m is obtained. The same reasoning, however, applies also to all other terms in 
the sum in m-- for each of the terms with |n| > 1, one can obtain its leading exponential 
contribution by performing exactly the same steps that Liischer did for the |n| = 1 term and 
work out three of the four integrations explicitly. It is easy to keep track of the vector n in 
doing these manipulations, and to get the resummed formula 

Mp{L) -Mp = - Y ^ r + 0(e-^^) (19) 

3271^ Xp y/n J-oo 

where m{n) has been given in Tab. Hand Tp{i') is the Pn forward scattering amplitude as usual. 
The extension which we are proposing is done in the same spirit as the extension of the domain 
of integration in (HI to inhnity (the contributions from the region ||/| > are beyond 

the accuracy of the formula): being of “kinematical” nature, the extension comes at no cost 
and may be numerically relevant. In this case, actually, one even obtains an improvement in 
the algebraic accuracy of the formula: we now have M > and not, as before, M =\/2Mt^. 

In other words, the terms and 0(e“^^’") are included now. To further clarify the 

meaning of (HnD let us list the two main classes of exponentially suppressed contributions which 
are still missing: 

1. All diagrams which have more than one pion loop in hnite volume. Obviously, these 
contributions start at the two-loop level in the chiral expansion. 

2. Contributions to the integral m which are due to singularities in either the propagator 
Go or the vertex function T which are further away from the real axis than These sin¬ 
gularities show up only if one considers the vertex function at one loop, or the propagator 
at two-loop accuracy and beyond. All these contributions appear only if one calculates 
hnite-volume effects at the two-loop level in the chiral expansion. 

In the hrst class we distinguish those diagrams where different loops factorize (i.e. loops which 
have no propagator in common) from those which do not. It is easy to see that the former 
sub-class yields corrections of order For two-loop diagrams which do not factorize (e.g. 

the two-loop sunset diagram) on may rely on the general discussion by Liischer (cf. in particular 
Eqn. (2.49) in Ref. [12]) and conclude that the sunset diagram decays exponentially at large L 
as exp(—Asunset^jr) with Asunset = {^3+1)/\/2 ~ 1.93. For the second class we remark that the 
singularities neglected in (HnD are due to the exchange of at least two pions, hence starting at 
z/ = ±2M.n-, and this gives terms of order We have not tried to prove the statement that 

the algebraic accuracy of our formula is given by 

M = M^{V3 + l)/V2 (20) 
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beyond the two-loop level, since it appears to us to be a question of academic interest. 

More interesting is the question of how accurate the formula (P|) is numerically. A complete 
analysis of hnite-volume effects for and iA at the two loop level, which is currently under 
way 1201, will clarify this point. The partial results we have so far indicate that the formula is 
very accurate. We have tried to hnd an algebraic reason for this, and found out that at the 
two loop level all diagrams which do not appear in Eq. (d are suppressed (besides an extra 
exponential factor) by some power of 1/L. The numerical results, however, seem to go beyond 
what one would expect from such an argument. 

An analogy to the low-temperature expansion seems more suggestive. In the effective theory 
the large volume and the low temperature expansions are in one-to-one correspondence |3]. In 
Ref. j2I] Schenk discussed the propagation of pions through matter in a state of thermal equi¬ 
librium at inverse temperature /?. If the temperature is not too high, the hadronic phase mainly 
consists of pions, with effects of other excitations such as K^rj^p... exponentially suppressed. 
Due to interactions with pions of the heat bath, the effective pion mass M-,^{j3) is given by j25 

Mm - M. = - 5 ^ nm) m + 0(nl) (21) 

with ujq = \J and the density n-Q^x) = —1). The details of the pion kinematics will 

be discussed in the next section, and the isospin index refers to the f-channel. One immediately 
verihes that m agrees with the modihed Liischer formula in one dimension to hrst order in the 
density ne - the density factor ub in the integrand is the outcome of the resummation over n in 
(fTR|l if n is taken as a one-dimensional vector. Schenk has carried the expansion of m one step 
further and determined the contributions of order to the pion mass at hnite temperature: in 
this extension, effects generated by three-body collisions are explicitly accounted for. It turned 
out that these effects are numerically very small, in line with intuition - the rescattering of 
three pions into three pions is a rare process unless the density is very high. Although the 
argument cannot be formulated in the same way for the hnite volume case, we do see that the 
outcome of the numerical analysis is the same. 


4 Meson masses and decay constants in finite volume 


We start with the resummed Liischer formulae for the relative hnite-size shift of pseudoscalar 
masses and decay constants 


p — 

RiPp — 


Mp{L) - Mp 




Mp 

Fp{L) - Fp 


Fp 


V r dy + 0{e-^^) (22) 

t^rdyAfp{iy)e-V^)^^ + Oie-^^) (23) 
IbMrpAp ^ ^/n J-oo 


327r2MpAp — 

m{n) 


where all symbols on the r.h.s. refer to inhnite volume quantities. The multiplicities m{n) are 
given in Tab.^ Ap has been dehned in m, and the dimensionless integration variable (which 
we will be using from here on) relates to the previous one via y = y/M^^. The amplitudes 
jFp(z/), A/p(z>) with i) = vlM.„ are 

Tp{y) = Tif{b),-^Mpv) 

Mp{i)) = -iA^=°(0,-4Mpz/) (24) 

where T^p^{t, u—s) is the PvT-scattering amplitude with zero t-channel isospin, and u—s) 

is the subtracted amplitude for the decay of the meson P with momentum p into two pions in 
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an isospin zero state (in the t-channel) via an axial cnrrent insertion. The snbtraction removes 
the one-particle redncible contribntion and is dehned throngh 


/ T/= 0 \ _ / tJ- 0 ] p '^tiF — s) 

)^l - Ih/M^P g2 _ jVf2 

= {{<PiXP2))i=M,mP{p)) (25) 

with Q = p — Pi — P 2 - Notice that the axial cnrrent in (j2ni) mnst be normalized snch that 
(0|y4^(0)|P(p)) = ipfiFp. The amplitndes Fp{9),Afp{u) for P = have all (with the 

exception of see below) been calcnlated at least at the one-loop level in ChPT. They 

have the generic form 

Pp(n) = FP{P) + ipF]l:\P) + epPP{^) + 0{ep) 

Mp{v) = Mp\v) + ^pN'p\v) + i^pN'p\v) + 0{^y) 


with defined in m- Inserting snch an expansion of the amplitnde in 


leads to 


Pm p 

Rpp 


_ ^ m{n) 1 Mtp 

2^/n KMp^^ 
_ ^ rn{n) I ^ 

~ h ^^Fp^- 


+ ^pi'mI+ epizi +o(«p) 

( 2 ) 


r(4) 


r2 AG) 


Fp’ + ipIK + 


(26) 

(27) 


where the I^FF) pg written in terms of a few basic integrals, as reported in Sect.El 

Please note the relative factor 2 in these two eqnations, to be consistent with da ESI In the 
following we elaborate on the explicit form of (dl for the pion, kaon and eta. 


4.1 Pion 

The amplitudes , J\f.p{A) defined in Eq. TIM have been given in Refs. daiia, respectively. 
We give them here for convenience. 

4.1.1 Pion mass 

Consider (Minkowski space) TTTr-scattering 

7r(pi) -h 7r(p2) ^ 7r(p3) -F 7r(p4) (28) 

with the kinematics 

s = (pi + , t={pi- p^f , u={pi- pif . (29) 

Isospin decomposition allows one to construct the f-channel isospin zero amplitude 

u- s) = 4(s, t, u) -1- 3A{t, s, u) + A{u, s, t) (30) 

from the invariant amplitude A{s, t, u) |221- amplitude entering the Liischer formula follows 
by imposing the forward scattering kinematics f = 0, viz. 

P^(z>)=T/=°(0,-4M^i/) . (31) 
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4.1.2 Pion decay constant 

We adopt the notation of Ref. usi, which relates to the one used in Eq. (I2S1) by means of the 
one pion in the initial state being transferred to an outgoing pion. Crossing symmetry relates 
the two via ps = —p. The amplitude for the creation of three pions out of the vacuum with an 
axial current has been performed up to NLO in Ref. [22] • It is decomposed according to 

(7r^(pi)7r^(p2)vr^(p3)|^il(0)|0) = {pi + P 2 )G + {pi - P 2 )H + p^^ F (32) 

with the three scalar functions ^ F = F(si, § 2 , ss), G = 82 , 53 ) and H = H{si, S 2 , S 3 ). 

The superscripts on the pion states and axial current are isospin indices. We have employed 
the variables Si = {p 2 + Ps)^ and cyclic permutations. The combination which has two of the 
outcoming pions in an J = 0 state in the S 3 -channel is given by 

((vr(pi)7r(p2))r=ovr=^(p3)|dJ(0)|0) = (df°), 

(d[=°)^ = (pi + P 2 )/, G'o(si, S 2 , S 3 ) + (Pl - P 2 )^, Ho{si, S 2 , S 3 ) + (P 3 )^ Fo(si, S 2 , S 3 ) (33) 

with the isospin projected ^ form factors 

To(si, S 2 , S 3 ) = 3Fi23 + G 23 I + G 312 — H 231 + H 312 

Go{si, S 2 , S 3 ) = 3Gi23 + 2 [-^231 + T 312 + G 231 + G 312 + H 231 — H 312 ] 

-^^0(515 S2, S3) = 3 Hi 23 + 2 [-^231 ~ -^312 ~ G231 + G312 — H231 — H312] 

where Xij^ = X{si, Sj, Sk) and X = F,G,H. The pole that the amplitude (A^"^)^ has in the 

unphysical region Q^ = M^, Q = — (pi +P 2 + P 3 ), needs to be subtracted as specihed in (OK|l 

(A'-"), = (A’-X - ~ (34) 

and in the end the result is evaluated in the forward kinematic conhguration, i.e. for S 3 = 0 . 
Hence the function is a function of just one variable i^ = (s 2 — si)/(4Mn-), viz. 

Hf°(0,-4M^z/) = 2vho{v) + Mj,{v) (35) 

^o(^) = hfo( 2 M 7 r(M^ — z/), 2 M.n-(M.n-+ I/), 0 ) 

/o(z/) = Fo{2M^{M^ - z/), 2Mt,{M^ + i^), 0) 

where the bar on the Fq form factor indicates that it is dehned after the subtraction of the pion 
pole (TTo remains unaffected). The amplitude which enters the formula for Fn-(T) is then 

A4(Z>) = -iHf°(0,-4M^z/) . (36) 


4.2 Kaon 

The amplitudes Fk{,^) and Nxi,^) niay be extracted from the form factors of the F ;4 decay, as 
calculated in Ref. [221 up to NLO. Furthermore, an explicit representation of Fi^(h) is found in 
the TriL-scattering study of Ref. |2Hj which has been extended to NNLO in Ref. [221 • Below, we 
stick to the NLO input, both for Fi^(z/) and A/i^(i 2 ). For the mass the hnite volume corrections 

^Note that our functions F, G, H are l/-\/2 times those of Ref. |22|, in agreement with the notation in Ref. m 
^The relation for Gq is given for completeness; on imposing forward scattering kinematics, it will drop out. 




at 2-loop level are very small (see discussion below) and a NNLO input is therefore not of 
particular interest for a lattice application. For the decay constant the result of the Liischer 
formula with NLO input is not particularly small (see below), and a NNLO rehnement would 
be useful. However, in this case one of the ingredients (to be precise: one of the 2-loop form 
factors of the axial-vector matrix element discussed in M) is missing. In the following we 
establish the relation of the amplitudes and to results given in the literature. 

4.2.1 Kaon mass 

The vriF-scattering amplitude T.^K{s,t,u) has been calculated at NLO in Ref. 

ir+(pi) + K*(p2) !r+(p3) + (37) 

with the Mandelstam variables (I2ni). This gives the t-channel isospin zero amplitude via 

T^K^(t,u- s) = ^[T^K{s,t,u) +T^K{u,t,s)] ( 38 ) 

from which the amplitude entering the Liischer formula follows by applying the forward scatte¬ 
ring kinematics, t = 0, viz. 

.F^(z>) = Ti^°( 0 ,- 4 Mkz/) . ( 39 ) 


4.2.2 Kaon decay constant 


We need the matrix element of an axial current between a kaon and two pions, which occurs in 
the evaluation of the Ki^ decay. Ref. [21] dehnes the three scalar form factors F, G, R through 


(Ak), = ^(vr+(p0vr-(p2)K‘^(0)|K+(p)) 


{Ak)i, = 


{pi + P2)fM F + {pi- P2)^, G + Q^R 


with F = F{s,t,u),G = G{s,t,u), R = R{s,t,u) and the kinematic variables 

s = {pi+p2f, t = {Pi-pf, u = {p2-pf, Q = p-Pi-P2 
The combination with the pions in an s-channel isospin zero state is 


(40) 


(41) 




7=03 
K )l^ 


V2M, 


K 


{pi + P2)tJ. -l- {pi — P2)fiG -|- Q^R'^ 


where 




(42) 


with X = F,G, R. Subtracting the pole at Q'^ = {p — Pi — ^ 2 )^ = as dehned in (OK|) and 
evaluating the amplitude — u) in the forward scattering conhguration s = 0 yields 


i^=°(0,-4M,,z/) = 


3i 


V2M, 


K 


2ug (u) + 


(43) 


'*As noted in the literature, there are two typos in (3.16) of [21]: The prefactor of should 

read and the term multiplying |Tk,,(n) is (n—|(M^-|-M|-))^ [the latter is correct in the preprint]. 
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with 


g-{v) = G-(0, Ml + Ml- 2 Mku, + 2Mki^) (44) 

and analogously for . Here, the bar on the form factor indicates again that it is dehned 

after subtraction of the kaon pole (the form factor G~ remains unaffected by the subtraction). 
Finally, according to Ea. (12411 . the amplitude entering the Liischer formula is 

Mk{^) = -iH^=°(0, -AMku) . (45) 


4.3 Eta 

The decay constant of the r] is of no phenomenological interest, like for the other neutral 
members of the pseudoscalar octet. We therefore refrain from discussing the hnite volume 
effects for this quantity (as a side remark, we notice that also the analogue of the decay 
amplitude, the ((27r)j-=o|^^lh) amplitude, is of no phenomenological interest and has never been 
calculated). We restrict ourselves to the hnite volume effects on the r] mass. 

The Trr^-scattering amplitude Tjr^(s,t,M) has been calculated to NLO in Ref. ra, 

’r°(pi) + '?(P 2 ) ^ ^“(ps) + v(Pi) (46) 

with kinematic variables (I2ni), and relates to the isospin zero amplitude in the f-channel through 

s) =3Rnis,t,u) . (47) 

The amplitude entering the Liischer formula follows by applying the forward scattering kine¬ 
matics, t = 0, viz. 

= p-"(o, . (48) 


5 Summary of the analytical results 

In order to use the asymptotic formulae we have to feed them with the specihc expressions for 
the scattering amplitude or the axial vector matrix element that are available in the literature. 
In this section we present such explicit formulae for Mp{L) and Fp{L) in two versions. We 
start with the complete expressions with some of the lengthier parts relegated to the appendix. 
The second step entails simplihed versions of the unhandy parts, together with a discussion of 
how they relate to the complete version. 


5.1 Full formulae 


Evaluating (1.41 D in the SU(2) framework the fractional shift of the pion mass takes the form 

m with 


7'(2) _ 

- 

r(4) _ 

— 

r(6) _ oC 

^ 


^8- 5 

18 3 2 


rrh - 2L 


10049 13 ,, 20- 

- N H- G 

1296 72 9 


40 


+ R2 

3^ 


27 ^ 4 ^ 


112 

L~9~ 

110 


32 






?! 


5l| 


oo _ _ yn oo _ 47 

+ (l6£i -l- — 11^3)f*4 + + ^2li + ~ -^3 + ^4 + 


3 . " V 9 

-|-5ri -I- 4 r 2 -|- Srs -|- 8 r 4 -|- + IGrg 
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+B^ 

+^7r( 


3476 

100 


77,, 32- 464- 448. 32,- -,- 

N + —ii -^^2 + “^^4- 3 ”^^^ ^ 4£2)^4 

248 

4 -- 8r3 - 56r4 - 48r5 -h 16r6 


8- 176 

+ 3^‘ + “*■= 


9 3 9 9 

where we use the abbreviations (77* denotes the modihed Bessel function) 


+ ^' 


( 6 ) 


(49) 


B\x) = 2Ki{x) , B‘^{x) = 2 K 2 {x)Ix 


ip = ln(M|.//x2), 


N = Wtt^ 


(60) 

(51) 

(52) 


and the ii which carry a mild logarithmic quark mass dependence m 

l. = ff>'' + 21og(-^) 


(53) 


with mass independent if^^^ as given in Tab. El The terms are contributions from 

the loop functions at order They are explicitly given in appendix^ and their numerical 

importance is discussed below. Formula (im|l has already appeared in jHj. 

The fractional shift in the pion decay constant takes the form (EH) with 


I^p'> = -25° 


4? = 


7-4- 
■- + 2£i + -4 
9 3 


374 


+ B^ 


112 

~9~ 


32 


+ S 


(4) 


(54) 


and S'^ moved to the appendix. Formula dSl has already been given in na. 

With the abbreviation xpq = Mp/Mq the hnite volume shift of the kaon mass is given by 


r(2) 

^Mk 

r(4) 

^Mk 


= 0 




ttK 


^ + 8 NxM^L\ + LI 


4L^ 


LI + 4L^ + 2 LI) 


+7j, 


x: 


ttK 


+ 




(- 


1 - lOx^ri^ + xIk 7 + X.,rK ' 


4(1 - X^k) 16 ^ 1 - X^K Q{Xr,K - 1) 


T or ) f Q (7 Xtjk){x^p[ 7) T r . Xt^x 


53 


X 


ttK 


32 V 3 ■ 5- /V-'/- -/ ■ g 

+52 -8Nx^k{^LI +LI)-i. 

Xq^K 


1 - lOx^K + xIj^ ' 

3{Xr,K - 1) 






1 ^TjK 1 


2(1 - x^k) 


^ 2 {x^k - 1). 


-I 0 ( 4 ) 


and the hnite volume correction for Fk takes the form (I2ID with 


(65) 


7'(2) _ 

^Fk - 

74 


--S° 

4 


16 


^ttK 1 


X-jyK^tt ^ ^r) XXr^^K ^ 2 ^ ^ 


X. 


XjtK 


XKt} 


Tvrj 


+ + 8xnKir]) + ^2NxTtK{,^L\ + 


2L\ 


-3)) 

3iVL5(l + x^i^) 


+ B'^x^k 


15 i]^ Xj^J^i^ ^ 3 irq Xj^rqip^ 

-2 1 X'j^x 2 1 Xxyi 


2AN{ALl + Ll) +S 


(4) 

Fk 


(56) 


(57) 
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and Sp^ given in App.|^and the convention (jKn|l applied thronghont, i.e. lx]. = -^x].(v^-^ 7 r) 
and Sxl, = Sxl,{y/nXT,). Analogonsly, the eta mass in hnite volume is given through 


7'(2) _ 

^Mr, - 




r(4) 


= xiLHB^ 


nr} 


2 + X. 


nrj 


+ t-^x 


( ---— ) 


n'^nrj 


'3(1 “^nr}) 
2 


+t'x(2a;x^ - - - 

^ 6 3(1 

+16Ar(6(L^i -Ll + Ll- L\) + + x^r,{<oL\ + 3L^ 

+ 5 ^ [9(1 + ix) - 32N{3Ll + LI) 


+ ^mI 


(58) 


Comparing dnH) and (P|) to (P|l it is obvious that the SU{3) breaking renders the expressions 
for Mk{L), Mj^{L) substantially more complicated than for It is remarkable that to 

leading order the hnite volume correction for the kaon mass vanishes [in the theory with virtual 
pions only, cf. Eq. (innD in app. P, while for the eta mass there is a suppression factor (M^/M^)^ 
[cf. Eq. flHH]) in app.p. As we shall see in the numerical analysis, these hnite volume corrections 
are practically negligible. 


5.2 Simplified formulae 

Beyond tree-level, the chiral representation of the amplitudes jFp(h),A/p(z/) that enter the 
resummed asymptotic formulae tend to become rather complicated. As a result, the expressions 
and Sp^ are not particularly handy, see App. |3 In the immediate vicinity of z/ = 0, 
however, a polynomial approximation to the chiral amplitudes reproduces them rather well. 
The reason behind is that the nonanalytic structure closest to the origin is the cut starting at 
z/ = ±Mjr. Therefore, for imaginary u close to the origin [i.e. for |/G [—M,r, Afn-] with y relating 
to V via m] even a second order polynomial reproduces the amplitude rather accurately, while 
outside this region the quality of the representation does not matter, due to the suppression 
factor exp(—Y^n(M2 -1- y‘^)L). The advantage of such a polynomial representation is that all 
integrals can be performed analytically, and there is an analytic bound on the remainder. 

For the pion mass and decay constant we hnd 

- ^(Af)go + 32g, + 2 ^ 32 ) 3 ^ + 0(b\^K)) 

= ^(8(70 - 13 ^ 1 ) 5 ° - ^(40^0 - 12(71 - 8(72 - 13^3)i?' + 0{b\^K)) (59) 

where the coefficients gi are the Taylor coefficients of the function ® 

^(x) = CT log + 2 (60) 

a -I- i 

with a = — 4/x around the point x = 2 

fi'(2 + e) = 5'o + + 7:92^^ + + 0{e^) (61) 

_ 2 6 

^The function g{x) relates to the standard scalar one-loop function J through g{x) = 16 Tr^J{xM^). 
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i 

L'(M,). 10’ 

1 

0.38 ±0.18 

2 

1.59 ±0.15 

3 

-2.91 ±0.32 

4 

0.00 ±0.80 

5 

1.46 ±0.10 

6 

0.00 ±0.30 

7 

-0.49 ±0.24 

8 

1.00 ±0.21 

9 

6.90 ±0.70 


i 

ri{Mp) 

1 

-1.5 X (1 ± 1) 

2 

3.2 X (1 ± 1) 

3 

-4.2 X (1 ± 1) 

4 

-2.5 X (1 ± 1) 

5 

3.8 ± 1.0 

6 

1.0 ±0.1 


i 

£phys 

1 

-0.4 ±0.6 

2 

4.3 ±0.1 

3 

2.9 ±2.4 

4 

4.4 ±0.2 


Table 2: S'f/(2)-framework: values at the physical pion mass of the low energy constants ^i 
from [21] together with the p® low energy constants fi{^ = Mp). Note that the £i used in the 
formulae for differ from these values by a term logarithmic in , see (El. 

S't/(3)-framework: L\{^ = Mp) taken from the 0{p^) fit in [20] (the uncertainties in Lg and Lg 
are our estimate). 


with the explicit values 


71 71 1 1 71 Svr 1 

50-2--, 5i-^-2> ~ 2~ 8 ' ~ ~W ~ 2 


(62) 


For 5^^ a similar short-hand version follows in the same manner. We refrain from showing 
them here, because these contributions turn out to be numerically so small that one could 
simply drop them. In and a polynomial expansion would not really simplify the 

representation, due to the different meson masses involved in the loop functions. 


6 Numerical analysis 

We are now in a position to evaluate the formulae for the relative hnite volume corrections to 
Mp and Fp, as presented in the previous section. To fully specify the meaning of our formulae 
we need to give, as the last ingredient, the quark mass dependence of the inhnite-volume 
quantities M^-, M^, iv, In line with the setup of our calculation we use 2-flavor ChPT 
for the pion mass and decay constant and 3-flavor ChPT for the kaon and eta counterparts. 
In either case the low energy parameters are determined from phenomenology and summarized 
in Tab. El Regarding the SU{2) low energy constants It and rj(p) we use the values obtained 
in Ref. 1231221; for the SU{3) low energy constants T-(h) 'w® fo the O(p^) fit of Ref. [HU] . 
In the former case, the full correlation matrix is given in j22] in the latter case it has been 
communicated privately |Sl|, but in either case the hnal errors are almost the same, regardless 
whether the full correlation matrix is used or just the diagonal part. 

6.1 Mtt dependence of 

The quark mass dependence of and F^^ has been computed, up to 2 loops, in Refs. j22lESll2n]- 
What we need in the present context is TV as a function of M^r, i.e. the single relationship that 
one gets after eliminating the quark mass. This relation has been given in dl, where also 
the SU{2) low energy constant F = (86.2 ± 0.5) MeV has been found. We stress that with 
Mjr (or iv) we mean simultaneously the pion mass (decay constant) in an infinite volume 
lattice simulation and in ChPT to the highest loop order available, i.e. to 0{p^) in the SU{2) 
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framework. Finally, we mention that the chiral expansion parameter = (^7r/4vrFV)^ remains 
small for pion masses up to 500 MeV, see the discussion in Hg for details. 


6.2 Mjr and dependence of Fk and 

In the 3-flavor case we have two independent quark masses, the average down and up quark 
mass and the strange quark mass. We take the liberty to rewrite everything in terms of 
and nis, for reasons that will become obvious soon. Then the 1-loop quark mass formulae read 


with 


= M| + 


= M'^ + 




M^{ - 2k 

+ 'msBo)^S{ki + 2 ^ 2 ) + 


F2 


+ 2 ^ 3 ) + 

o c 64 4 o 2 o 'i 

+M^msBo^ — {ki + 3^2 — 2ks) + in — - ^r?)} 

, T-> \2( ^28, 3 , 4 o 8o'|] 

+ {msBo) {-^(^1 + 2^2 + h) + 


Fx — Fh -l- 


4(Mi - M^L' + - jMIIk - ^MX} 


AW8 


k, = 2Ll - LI , k2 = 2Ll -LI, k3 = 3Uj + U, 
and with the abbreviations 

Ml = m.Ba + ^Ml , = l(Mj + 4m.B„) 


(63) 


(64) 

(65) 


( 66 ) 


together with N and ip as dehned in (1^ and accordingly ip = ln(M|,/p^). Note that M|- 
and are of hybrid nature ~ the part refers to 1-loop ChPT, while the part linear in 

rris is a tree-level contribution. This is unavoidable if we want to discuss the dependence of 
physical quantities on the pion mass, instead of on quark masses. In technical terms (insD is 
used to £x rrisBQ] we simply require that Mx takes the physical value for the 

result is BquIs = 0.223GeV^. The analogous requirement for implies that we must slightly 
readjust L\ to —0.47 ■ 10“^ [well compatible with the error in Tab.|2], for all other low energy 
constants the central values in Tab.|2lare used. Even in the 3-flavor case we choose to describe 
the quark mass dependence of and through SLf (2) ChPT, as discussed in the previous 
subsection. Note that for this choice exactly reproduces what one would get in the 

SU{3) framework, since the phenomenological ii values know about the virtual strange quark 
loops. In actual lattice simulations rUg is typically close to the physical value, and we expect 
this to remain a valid approximation. The resulting ikC dependence of Fx, Mx, is shown 
in Fig.^ One notices that Mjj ~ 640 MeV for ikC ~ 500 MeV, thus the expansion parameter 
^P remains small in the entire mass range < 500 MeV, even for P = r]. In our numerical 
analysis we will use ^p exactly as determined from Fig.^ and ignore the uncertainty of this 
computed expansion parameter, since in a lattice computation one may iteratively determine 
Mp{L) and Ft^{L) and thus ^p. We do, however, consider the uncertainties in the expansion 
coefficients Ip^f* in and (j27j), respectively, with details specihed below. We 

have checked that even including the contribution of ^p to the total error would barely change 
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Figure 1: dependence (in infinite volume) of Fk, Mk, For the latter three quantities 

the strange quark mass has been fixed as to reproduce at M^ = . 


the errors in our main result, Figs. l2ll5l In summary, the quark mass dependence of a quantity 
X = Fki Mk, Mn is considered a function of alone through 

X = X(M^, rUsBo = 0.223 GeV^) (67) 

and an appropriate choice of the SU{3) low energy constants, as given in Tab. El 

6.3 Results 

We plot our results for Rm^ in Fig. El both for L = 2,3,4fm as a function of (top) and for 
Mj, = 100, 300, 500 MeV as a function of L (bottom). The result of the original (“u = 1”) Liischer 
formula m with LO/NLO/NNLO chiral input is shown with a thin dotted/dashed/full line, 
respectively. The resummed (“all n”) formula (O with LO/NLO/NNLO chiral input is given 
with a thick dotted/dashed/full line, respectively. With NLO or NNLO input the pertinent low 
energy constants lead to a non-negligible error band, except for where the error with NLO 
input is of the order of the thickness of the line, with LO input it is zero. This is a consequence 
of our choice to disregard any uncertainty in the expansion parameter (here ^.n-), as discussed in 
the previous subsection. The area which corresponds - in the resummed scenario with NNLO 
input - to a situation with M.,^L < 2 should be disregarded, since there is no reason to hope 
that the resummed Liischer formula would still capture the numerically dominating terms in a 
complete “ChPT in finite volume” formula to the corresponding order in the p-counting. 

Fig. El contains our data for —Rf^,- They are organized in the same manner as in the previous 
figure, though only LO and NLO input is used, since the pertinent matrix element is known 
only to 1-loop order, as discussed in Sect. 4. 

Finally, Figs. lHIKl contain the same information for Fk, M^k, Fl^q. In all these cases the 
asymptotic formulae with and without resummation may be compared, and the effect of going 
from LO to NLO input may be assessed. 

The main message to be extracted from Figs. EE is that the relative finite volume shift 
vanishes indeed in proportion to in the logarithmic representation one has an almost 

linear fall-off pattern, and higher orders mainly affect the prefactor. The second point concerns 
the relative importance of higher orders in the chiral counting versus higher exponentials. For 
small pion masses (say 100 MeV) resumming (i.e. higher exponentials) prove vital, while for 
large pion masses (say 500 MeV) higher loop orders prove more relevant (though the effect is 
small in that regime). As has been discussed in Ref. ^1] a large shift in Rmp, when moving 
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(GeV) 



Figure 2: Rm^ vs. for L = 2,3,4fm (top) and vs. L for = 100, 300, 500 MeV (bottom). 
The result of the original (“n = 1”) Liischer formula with LO/NLO/NNLO chiral input 
is to be compared to the resummed (“all n”) formula (ITl) which amounts to an approximate 
1/2/3-loop ChPT calculation in hnite volume. With NNLO input the low energy constants 
lead to a non-negligible error band; with NLO input the error is smaller (not shown), with LO 
input it is zero (see text). In the region above the Mt^L = 2 line one is not safely in the p-regime 
and our results should not be trusted. 
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Figure 3: —Rf^ vs. for L = 2,3,4fm (left) and vs. L for = 100, 300, 500 MeV (right). 



Figure 4: 


• • • • LO, n = 1 

- - NLO. n = 1 
■ - LO, all n 

— NLO, all n 
_ ML = 2 




M,j(GeV) L (fm) 

-Rfk vs. Mjt for L = 2,3,4fm (left) and vs. L for = 100, 300, 500 MeV (right). 




Figure 5: Rmk (left) and Rm,, (right) vs. M^- for L = 2,3,4fm. 
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Rm^ 
140 MeV 
160 MeV 
180 MeV 
200 MeV 
220 MeV 
240 MeV 
260 MeV 
280 MeV 
300 MeV 
320 MeV 
340 MeV 
360 MeV 
380 MeV 
400 MeV 
420 MeV 
440 MeV 
460 MeV 
480 MeV 
500 MeV 


1.6 fm 
.2099(127) 
.1687(130) 
.1366(130) 
.1113(127) 
.0912(123) 
.0749(117) 
.0618(110) 
.0511(102) 
0.0423(93) 
0.0352(85) 
0.0293(77) 
0.0245(69) 
0.0205(62) 
0.0172(55) 
0.0145(49) 
0.0123(43) 
0.0104(38) 
0.0088(33) 
0.0076(29) 


2.0 fm 
0.0830(55) 
0.0647(55) 
0.0509(53) 
0.0403(50) 
0.0320(47) 
0.0256(43) 
0.0205(39) 
0.0164(35) 
0.0132(31) 
0.0107(28) 
0.0086(24) 
0.0070(21) 
0.0057(18) 
0.0046(15) 
0.0038(13) 
0.0031(11) 
0.0025(10) 
0.0021(08) 
0.0017(07) 


2.2 fm 
0.0552(38) 
0.0424(37) 
0.0328(35) 
0.0256(33) 
0.0200(31) 
0.0157(28) 
0.0124(25) 
0.0098(22) 
0.0078(19) 
0.0062(16) 
0.0049(14) 
0.0039(12) 
0.0031(10) 
0.0025(09) 
0.0020(07) 
0.0016(06) 
0.0013(05) 
0.0011(04) 
0.0009(03) 


2.4 fm 
0.0377(27) 
0.0285(26) 
0.0217(24) 
0.0167(22) 
0.0129(20) 
0.0099(18) 
0.0077(16) 
0.0060(14) 
0.0047(12) 
0.0037(10) 
0.0029(08) 
0.0023(07) 
0.0018(06) 
0.0014(05) 
0.0011(04) 
0.0009(03) 
0.0007(03) 
0.0006(02) 
0.0004(02) 


2.6 fm 
0.0264(19) 
0.0196(18) 
0.0147(17) 
0.0111(15) 
0.0084(14) 
0.0064(12) 
0.0049(10) 
0.0038(09) 
0.0029(07) 
0.0022(06) 
0.0017(05) 
0.0013(04) 
0.0010(03) 
0.0008(03) 
0.0006(02) 
0.0005(02) 
0.0004(01) 
0.0003(01) 
0 . 0002 ( 01 ) 


3.0 fm 
0.0136(10) 
0.0098(10) 
0.0071(09) 
0.0052(08) 
0.0038(06) 
0.0028(05) 
0.0021(05) 
0.0015(04) 
0.0011(03) 
0.0009(02) 
0.0006(02) 
0.0005(02) 
0.0004(01) 
0.0003(01) 
0 . 0002 ( 01 ) 
0 . 0001 ( 01 ) 
0 . 0001 ( 00 ) 
0 . 0001 ( 00 ) 
0 . 0001 ( 00 ) 


2.8 fm 
0.0188(14) 
0.0138(13) 
0 . 0102 ( 12 ) 
0.0076(11) 
0.0056(09) 
0.0042(08) 
0.0032(07) 
0.0024(06) 
0.0018(05) 
0.0014(04) 
0.0010(03) 
0.0008(03) 
0.0006(02) 
0.0005(02) 
0.0003(01) 
0.0003(01) 
0 . 0002 ( 01 ) 
0 . 0002 ( 01 ) 
0 . 0001 ( 00 ) 


1.8 fm 
0.1292(82) 
0.1023(83) 
0.0817(82) 
0.0656(79) 
0.0529(75) 
0.0429(70) 
0.0349(65) 
0.0284(59) 
0.0232(53) 
0.0190(48) 
0.0156(42) 
0.0129(37) 
0.0106(33) 
0.0088(29) 
0.0073(25) 
0.0061(22) 
0.0051(19) 
0.0042(16) 
0.0036(14) 


Table 3: via the resummed Liischer formula ()22|1 with NNLO chiral input for jF^(z/), 

representing an approximate 3-loop result. The error includes the uncertainty of the ^i and the 
0(p®) low energy constants, but no systematics. Entries with < 2 are unlikely to really 
capture the physical hnite size effect, and the hrst two columns are somewhat on the short side 
with respect to the condition 0 - 


~Rft, 


140 

160 

180 

200 

220 

240 

260 

280 

300 

320 

340 

360 

380 

400 

420 

440 

460 

480 

500 


MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 

MeV 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


1.6 fm 
.5844(43) 
.4551(36) 
.3580(30) 
.2840(25) 
.2267(22) 
.1820(19) 
.1467(16) 
.1188(14) 
.0965(13) 
.0786(11) 
.0642(10) 
.0527(09) 
.0433(08) 
.0356(07) 
.0294(06) 
.0244(05) 
.0202(05) 
.0169(04) 
.0141(04) 


1.8 fm 
0.3683(24) 
0.2828(20) 
0.2194(17) 
0.1715(14) 
0.1350(12) 
0.1068(11) 
0.0848(09) 
0.0677(08) 
0.0541(07) 
0.0434(06) 
0.0350(05) 
0.0282(05) 
0.0228(04) 
0.0185(04) 
0.0150(03) 
0.0122(03) 
0 . 0100 ( 02 ) 
0.0082(02) 
0.0067(02) 


2.0 fm 
0.2414(14) 
0.1827(12) 
0.1397(10) 
0.1076(08) 
0.0834(07) 
0.0650(06) 
0.0509(05) 
0.0399(05) 
0.0315(04) 
0.0249(04) 
0.0197(03) 
0.0156(03) 
0.0124(02) 
0.0099(02) 
0.0079(02) 
0.0063(01) 
0.0051(01) 
0.0041(01) 
0.0033(01) 


2.2 fm 
0.1633(09) 
0.1218(07) 
0.0917(06) 
0.0696(05) 
0.0531(05) 
0.0408(04) 
0.0314(03) 
0.0243(03) 
0.0188(03) 
0.0146(02) 
0.0114(02) 
0.0089(02) 
0.0070(01) 
0.0055(01) 
0.0043(01) 
0.0034(01) 
0.0027(01) 
0 . 0021 ( 01 ) 
0.0017(00) 


2.4 fm 
0.1134(06) 
0.0833(05) 
0.0618(04) 
0.0461(03) 
0.0347(03) 
0.0262(03) 
0.0198(02) 
0.0151(02) 
0.0115(02) 
0.0088(01) 
0.0067(01) 
0.0052(01) 
0.0040(01) 
0.0031(01) 
0.0024(01) 
0.0018(00) 
0.0014(00) 
0 . 0011 ( 00 ) 
0.0009(00) 


2.6 fm 
0.0804(04) 
0.0581(03) 
0.0424(03) 
0.0312(02) 
0.0231(02) 
0.0171(02) 
0.0128(01) 
0.0095(01) 
0.0072(01) 
0.0054(01) 
0.0040(01) 
0.0031(01) 
0.0023(00) 
0.0017(00) 
0.0013(00) 
0 . 0010 ( 00 ) 
0.0008(00) 
0.0006(00) 
0.0004(00) 


2.8 fm 
0.0581(03) 
0.0413(02) 
0.0297(02) 
0.0215(02) 
0.0156(01) 
0.0114(01) 
0.0084(01) 
0.0061(01) 
0.0045(01) 
0.0033(01) 
0.0025(00) 
0.0018(00) 
0.0014(00) 
0 . 0010 ( 00 ) 
0.0008(00) 
0.0006(00) 
0.0004(00) 
0.0003(00) 
0 . 0002 ( 00 ) 


3.0 fm 
0.0426(02) 
0.0298(02) 
0 . 0211 ( 01 ) 
0.0150(01) 
0.0107(01) 
0.0077(01) 
0.0056(01) 
0.0040(01) 
0.0029(00) 
0 . 0021 ( 00 ) 
0.0015(00) 
0 . 0011 ( 00 ) 
0.0008(00) 
0.0006(00) 
0.0004(00) 
0.0003(00) 
0 . 0002 ( 00 ) 
0 . 0002 ( 00 ) 
0 . 0001 ( 00 ) 


Table 4: —Rft, via the resummed Liischer formula (0!^ with NLO chiral input for A/' 7 r(h), 
representing an approximate 2-loop result. The error includes the uncertainty of the £*, but no 
systematics. Entries with Mt^L<2 are unlikely to really capture the physical hnite size effect, 
and the hrst two columns are somewhat on the short side with respect to the condition 0 . 
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from LO to NLO input, does not necessarily signal a bad chiral convergence behavior, since 
the cut in the underlying jFp amplitude starts only at the NLO level. Unfortunately, the effect 
due to an upgrade to NNLO input can only be checked in the case of Rm^i since only there 
the pertinent amplitude is known. Nonetheless, we believe that the regime in the 
plane that leads to a nice convergence behavior in Rmt, is indicative of the regime where the 
resummed formulae for Rp^, Rmk^ Rfki with NLO input yield a trustworthy result. 

For Rm^ and Rp^ the numerical results to the highest loop order available (equivalent to 
an approximate 3-loop and 2-loop calculation in ChPT) have been collected in Tabs.OJandEl 
respectively. From the general discussion it is clear that the logarithms of the numbers in these 
tables may be interpolated with a low-order polynomial. 

7 Two types of applications 

We hnish with a discussion of two prototype applications of our formulae. The hrst one is a 
“forward-type” application, in which our formulae are used to control a systematic error in a 
lattice calculation. The second one concerns a “backward-type” application, where one tries to 
determine QCD low energy constants from explicitly measuring hnite-volume effects. 

7.1 Finite volume effects and Marciano’s determination of 14s 

An example of how an analytic hnite volume calculation may help to control a systematic error 
is the following. Marciano pointed out that, modulo radiative corrections, the ratio is 

hxed by the ratio of branching ratios for K 12 and 7^12 decays. Taking into account radiative 
corrections, he obtained the following relation 03 

= 0.07602(23)(27) (68) 

\yud\ ^TT 

where the errors represent the experimental and radiative correction uncertainties. He then 
suggested to combine the value for Vud obtained from superallowed nuclear beta decays with 
a value for Fk/F-w from lattice simulations. We stress that the necessary accuracy to make 
an impact on the determination of Ws is at the level of 1% or better, and indeed, both the 
determination of Vud as well as the ratio of branching ratios are known to well below 1%. This 
means that any improvement in the lattice calculation of Fk/F-j^ will be immediately reflected 
in the value of Ws- In particular, being able to control systematic effects to well below 1% is 
of crucial importance. 

With our results for Rp^ and Rpj^ it is straightforward to calculate the hnite-volume shift 
of the latter ratio 

= ^ {l + flFK-flF, + 0(flP} (69) 

and thus to compare the magnitude of this effect to the typical size of the statistical error. A 
plot of the hnite volume ehect for the ratio of decay constants as a function of the pion mass 
and for a few volume sizes is provided in Fig. IHl 

In his analysis Marciano uses the MILC Collaboration result = 1.201(8)(15) [HKl 13fi] . 

Among the various data sets they have, those with the smallest M^-L and thus most likely to be 
ahected by sizeable hnite volume corrections, have Mn-(L) ~311 MeV, L~2.4fm and Mn-(L)~ 
262 MeV, L~2.89fm. Using (IHn]) we hnd that with the parameters of the hrst set Fk{L)/Ft,{L) 
in the continuum deviates from the inhnite volume result by 0.0099 — 0.0038 = 0.61%, while 
the corresponding estimate for the second set reads 0.0067 — 0.0025 = 0.42%. A systematic 
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(GeV) 


Figure 6: The relative finite volume effect {Fk{L)F.^)/{F.^{L)Fk) — 1 vs. for L = 2,3,4fm. 
Above the M-j^L = 2 line one is not safely in the p-regime and our results should not be trusted. 

effect of half a percent needs to be taken into account in a high precision study, and from Fig.lHl 
one sees that other (M^, L) pairs to reach that level would be (410 MeV, 2 fm), (230 MeV, 3 fm) 
and (150MeV,4fm). 

We stress that the numerical example just discussed is for illustrative purposes only, because 
we do not know whether our formulae can be applied to the MILC Collaboration data. Lattice 
QCD with staggered fermions and A^=2orA^=2 + l violates flavour (or taste) symmetry 
and low energy unitarity, properties our analysis relies on. Under the assumption that these 
effects disappear in the continuum limit the finite volume shift of an observable like Fk/FV 
can be calculated in staggered chiral perturbation theory (this is what the MILC Collaboration 
does [SHI), cannot enjoy the benefits of the Liischer formula. We stress that no such 

conceptual issues arise with dynamical Wilson-type fermions. In such a case our continuum 
formulae can be directly applied to the data at finite lattice spacing with cut-off effects bringing 
only mild (i.e. numerically irrelevant) modifications as discussed in App.lRl 

7.2 Low energy constants from finite volume effects 

In the present section we discuss whether one can use these finite volume effects to obtain 
information on the low energy constants from lattice calculations. At first sight this seems 
unpractical, because these effects are quite small and decay exponentially with M^F. This 
means that, roughly speaking, if one wishes to obtain information on (a combination of) low 
energy constants to a certain accuracy, one has to calculate the corresponding particle mass 
or decay constant to an accuracy which is about two orders of magnitude higher, and this is 
a challenge. The asymptotic formulae provide a connection between finite volume effects on 

®As of now, there is no proof, but rather a lively debate on this issue in the literature ETj. 
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two-point functions and (infinite volume) four-point functions, and thus give us access to low 
energy constants that appear as local contributions to four-point functions. The question is 
then what alternative ways one has to measure these constants on the lattice. It is well known 
that the Pvr scattering amplitudes which govern the hnite volume corrections to the mass of 
the P particle can be obtained more directly by evaluating the hnite volume dependence of the 
energy of the state of two particles P and tt enclosed in a box [23 • This method is more direct 
because the effect is suppressed only by powers of the volume rather than exponentially. Still, 
such a calculation is very difficult, also because the typical volume needed is quite large. 

There is another important difference between the two methods. As we have seen in Sect. 15.21 
the analytic representation of the hnite volume ehects is enormously simplihed if one Taylor 
expands the amplitudes in the Liischer-type formulae. Here, two terms in the Taylor expansion 
are enough for a very accurate representation. The extraction of the low energy constants from 
these ehects can be viewed as a two-step process: one hrst determines the values of the hrst two 
Taylor coefficients of the amplitudes at u = 0, and then from these the low energy constants. 
The chiral representation enters only in this second step. Analogously, if one uses the method 
with two particles in a box, one hrst determines the scattering lengths, and then extracts from 
these the relevant low energy constants. As discussed in j^H], the scattering lengths have a 
badly converging chiral expansion, such that only at very small quark masses one would be 
able to reliably extract the low energy constants. This happens because one is evaluating the 
amplitude on top of the threshold singularity. At z/ = 0, below threshold and away from any 
other singularity, the amplitude displays a better convergence. For all these reasons we believe 
it is worthwhile to explore this alternative route. 

The quantities which are worth considering for our scope are M^r, iA and Fk- We write the 
relative hnite volume shifts in the form 

Rx = R% + cx + f3‘^L\) ( 70 ) 

where represents contributions independent of the low energy constants, and the coefficients 
Cx are dehned as 

cm^ = , cf^ = , CFji = ■ ( 71 ) 

The functions are series of Bessel functions and depend on Att only 

^ ^ . (72) 

n=l V 

They are plotted in Fig. 0 The are the combinations of low energy constants that appear 
in each of the quantities to next-to-leading order. In particular we read off from the formulae 

(I22ES} and iHElEZD 



= 2(0 + fl2)-|4-4 

= 

-|(^“i + 40) 


= 2(0 + fl2)-3l4 

L% = 

-|(^“i + 40) 


= x^xi^Ll + LI- 2L\) - 1(1 + x^k)LI 

Lh = 

-2x^k{^LI + LI) 


It is interesting that both and FA are sensitive to the same combination of and ^2 (note 
that £3 and £4 can be pinned down directly from the quark mass dependence of and 
respectively). This means that one can determine the low energy constants which £x the tttt 
scattering amplitude also from the finite volume effects for FA where the tttt scattering amplitude 
does not appear. One should remark that the sensitivity to the combination (£1 - 1 - ^£ 2 ) is the 
same, whereas FA is a factor two more sensitive to the combination {£i + 4 ^ 2 )• 
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To discuss the numerics we fix = 300 MeV and L = 2 fm. In this setting j3^ = 0.22, 
/3^ = 0.10, and cm^ = —0.0026, = 0.0026. A change of one unit in the combination (£i + ^£ 2 ) 

then generates a shift of about 0.12% in both Mjr(L) and Ftj{L). According to Tab.Elthis linear 
combination is known from phenomenology to be 

(«1 + = -0-1 ± O-S . (73) 

On the other hand a change of one unit in (£1 + M 2 ) modifies F-„{L) by 0.07% and by 

half that much. This linear combination, however, is known to be larger from phenomenology 

(fi + 44)|„,.3„„„=9,l±0.5 . (74) 

This numerical example indicates that in order to get a reasonable account on these combina¬ 
tions of low energy constants, one would have to control the pion mass and decay constant for 
Mjr = 300MeV and L = 2fm to less than 1 permille, which is a real challenge. 

An alternative way, as already mentioned, would be to calculate the vrvr A-wave 1 = 2 
scattering lengths via Liischer’s method Isa There are some results with dynamical fermions 
for this quantity m, but not yet precise enough and for low enough pion masses that would 
allow an extraction of the relevant low energy constant. We mention in passing that both 
Oq and Uq (the latter is even more difficult to calculate on the lattice, and shows a very badly 
converging chiral expansion [22]) are sensitive to the combination £i + 2 (!. 2 i and therefore provide 
complementary information to the one that would be obtained from the finite volume effects 
we have discussed here. 

The numerics for the case is similar. Again, in order to get a sensitivity comparable to 
the one which characterizes the phenomenological determination (notice that the L\ constants 
are usually given in units of 10“^) one would have to calculate to the permille accuracy. 
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8 Conclusions 


In this paper we have carefully analyzed the hnite volume effects on masses and decay constants 
of the lightest pseudoscalar mesons. The theoretical framework in which we carry out this 
analysis has been set up long ago by Gasser and Leutwyler |1] and by Liischer [T2]. The 
former two discussed how ChPT can be adapted to the case of a hnite box and then be used to 
calculate hnite volume ehects. The latter author derived a formula valid for large volumes which 
expresses the hnite volume ehects in terms of an integral over a physical scattering amplitude. 
This formula does not rely on ChPT, but for QCD is best used in combination with ChPT: this 
is the tool to provide a representation of the necessary scattering amplitude as a function of 
quark masses. More recently a formula d la Liischer for decay constants has been derived by two 
of us HSj. In the present paper we have proposed a resummation of the asymptotic formulae as 
the best tool to study hnite volume ehects for masses and decay constants. The resummation 
is a simple, “kinematical” extension of the asymptotic formulae, which however does improve 
the algebraic accuracy of the formula. Used in combination with the chiral representation for 
the scattering amplitude in the integral it yields an accurate determination of hnite volume 
ehects. 

Our numerical results show that decay constants and masses of the pseudoscalars are in 
general little ahected by the hnite spatial size of the box (as soon as the box is large enough 
that ChPT can be applied, L > 2fm). For the smallest acceptable values of (in the p- 
regime, in which we are working, this quantity has to be larger than one) they are typically of 
the order of a few percent for M.^-, and Fk- We have seen that in the p-regime Mk and Mrj 
are practically insensitive to the box size. Independently of the exact size of these corrections 
we could always check the convergence of the chiral expansion and conclude that these hnite 
volume effects are under good theoretical control. This means that if one’s goal is to calculate 
masses or decay constants one can use the results of this paper to choose the volume in order 
to minimize the calculational costs. For example, by using a box of 2 fm size (which our results 
show to be sufficiently large) and explicitly correcting for the hnite volume ehects one can save 
computational costs with respect to a 2.5 fm size box which (for pion masses of about 300 MeV 
or larger) gives hnite volume ehects below 1%. The gain in CPU time that comes from such a 
reduction of the volume by almost a factor 2 can be used more fruitfully by pushing towards 
smaller lattice spacings and/or lighter quark masses. 

Since our results are theory predictions, an explicit check by lattice calculations would of 
course be very welcome, but it would require a high precision. Once this level of accuracy will be 
reached, our formulae will become particularly useful. First, in order to correct for these ehects 
in all applications which require a high precision, like the evaluation of the /F/ ratio which, 
as suggested by Marciano [31], leads to a determination of I4s jSHl- Second, if one wants to use 
these hnite size ehects as a means to determine low energy constants on the lattice. As we have 
discussed, despite the fact that these ehects are exponentially suppressed and numerically quite 
small in all practical situations, they do oher the advantage of involving in two-point functions 
low energy constants which, otherwise, appear only in four-point functions. The latter are 
quantities which are very difficult to determine directly on the lattice, and it may therefore 
turn out to be easier to see them indirectly, through a small correction to an “easy” quantity 
like the pion mass or decay constant. 
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A The integrals and S^pp 


In this appendix we give explicitly the contributions from loop functions to the hnite volume 
effects. These have been introduced and dehned in Sect. El The expression for the pion related 
integrals have already been given in mmsi. We give them here for convenience: 



where the integrals are dehned as 

C’ ^ jyyy'‘e-'/^^)>~g('\2 + 2\y) for (76) 

= I r dy f e-V”(i+S")^- iV^ K"(2 + 219 ) for I ^ (77) 

L im J—00 y k, ocicl 

with dehned in (Eni and 


K^ix) 

Krix) 

K^ix) 


1 r n 2 

g{x) - 2 


L- 

-^[g{xf - 4g{x) 
1 


2g{xy — 12g{xy + 24:g{x) + 2TT^a^g{x) — 16 — n^a^x 


2N‘^a^x f 

^[A'r(x) + iAT'(x) + iA7*(x) + 


(78) 


with a = — 4/x. As seen in fl76j) . in integrals over the derivative of the function g appear. 

This is a consequence of the subtraction procedure (13411 . In the practical implementation one 
writes si = 2M^ — 2vM.„ + — M ^)/2 — s^/2 and S 2 = 2M^ + 2 i/M^ + — M ^)/2 — S 3/2 to 

trade si, S 2 for z/ and expands consistently in these new variables. For instance 

in Fq one substitutes J(si) ^ J{2M^ — 2uMy + J'{2Ml — 2vMy[Q‘^ — Ml — S'^/2. 

The loop integrals for kaon and eta hnite volume ehects read: 


q(4) 

^Mk 


!{4(1 + x,a-)^SK - g(l + x,K)S]i. - fsil - 4(1 - 


X. 
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^(1 - - I + Sil) + 4(1 + 


96 

1 3 1 

-(1 + X^K)SjjK - g'S'^K + Y^(l + X^K)(3Xr,K + ‘^X^k “ 5)5'' 


0,3 

rjK 


^(1 - 2x^k + Xr,K)Sl^K 
Mk 


q(4) _ 


( cO.5 I c2,5\ 
[XttKO^X + ‘^tjK ) 


- 4(1 + ^^K)S\il - 


AU 
15 
16* 

^S’^Kxr + 4(1 - - 4(1 - XnK)S^’^ 


(1 + x^K)S]il + + ^(1 - 5x^k)S^/^ + (6 - ^x^k) 5 


1,3 

Ktt 


32 


16 


Ktt 


[ eO.5 I oc2,5 cO .6 c 2 , 6 \ 

[Xt.K^K'k + ~ X-kK^Kit ~ ^Kiz) 


9 
4 

3 

— (1 + XjtX ) 


1,1 


9 <^2,1 


— (1 + X^K^S^’^ + —(1 + X^k)S^k 


1,2 


9 


+ ( ^{XrjK - 1 ) - -^(1 + X^k) ) 5: 


16 


27 


32 


+ ^(4 - ‘2x^k + 3xr,K)Si’K + + ^(5 - 3a;^;^)(l - xlj^)S:, 


^0,3 


32 


0,4 

■’■qK 


2,3 

^ (1 - 2x^k + x,k)S^^’^ - I {x^KS^i + 2Sl^^ - x^KS^i - S^i) 


q(4) 

‘^Mr, 


16 


2 rpU,l rpl,l qrp. 

X-n-q-'-KK ’^X-j^q-LKK 


0,1 


+ -X T 

' ^‘^Trq qiT 


where we have introduced the following abbreviations 


XPQ — 


Ml 

Mi 


= In 


. F, 


(79) 


(80) 


The integrals are proportional to modihed Bessel functions 


iy = £(^(^)‘/4,,(v;7,,, (81, 


and the quantities S'pg and TpQ are integrals over functions qIq which occur at one-loop order 


nk,I 


in the chiral expansion. They are all analytical along the integration line. The expressions S'pg 


and TpQ are dehned as 


qk,I 

^PQ 


rjik,I 

^PQ 


with 


r dy f g%{Ml + Ml + 2 iMKM^y) 

im J—oo 


]\frik+l )/2 

Im 


fcp-V^hl^A. (/) 7 J^2 


y e 


gl>^{Ml + Ml+ 2 i M.My^y) 


for 


for ■ 


k even 
k odd 
k even 
k odd 


dpQix) = Jpq{x) , 9pQix) = MIJ'pq{x) 

fi'pQ(^) = ^pq{x) 5 9pq{x) = MIK'pq{x) 

9pq{x) = Mpq{x) , 9pq{x) = MIM'pq{x) . 

For completeness, we give the explicit expressions for the glq [21]. All functions can be ex¬ 
pressed in terms of the subtracted scalar integral J{t) = J(t)—J(0) evaluated in four dimensions 

d‘^p 1 


Jit) = -i/ 


(27r)'^ {{p + ky — M^){p‘^ — mX 


(82) 
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with t = The functions used in the text are then 


J{t) = 


J\t) 

m 

K\t) 

M{t) 

M\t) 

where 


^ dx In^ 


N Jo 
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tx{l 


X 
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A 2 
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- J\t) 


1 


1 - A 2 _ 

— {t- 2E| Jit) H- -Jit) H- 

12t ^ ^ ^ ^ 3t2 ^ ^ ^ 18A 


6Nt 


M'^mJ , w? 

^ + '"Mi 


1 
6^2 


^ J(«) + i(«^ - 2«E + 4A2) J'(i) + ^ (e + 


2M^m? mJ ' 


(83) 


A = M^ — rrr 


'L = M^ + nf 


p = p{t, M^, m^) = (t + A )2 — 4tM2 


(84) 


In the text these are used with subscripts 

JpQ{t) = Jit) with M = Mp,m = Mq (85) 

and similarly for the other symbols. We add a remark concerning the analyticity properties 
of the loop functions. The asymptotic formula requires them to be evaluated for complex 
arguments. There is one case, where the representation of Eq. P|) does not yet provide an 
unambiguous analytic continuation, namely for Jpqi^Mp + Mq + 2iMpMQy), because p = 
—4MpMq( 1 + p2). The correct analytical continuation is given by 

= 2iMpMquj ( 86 ) 

with 00 = a/ 1 + p2^ implying for the logarithm in Ea. (IHH|) (for t = Mp + Mq + 2iMpMQy), 

1 + y/P? - _ 1 ^ + y j +i7r for p < 0 

- v/p)'- A2 ’^cn-pj-ivT forp>0 ' 

All loop functions are now well dehned along the integration line in the asymptotic formulae. 


B Cut-off effects 

In this appendix we wish to discuss whether it is sufficient to calculate hnite volume effects in 
continuum ChPT or whether cut-off effects should be taken care of when correcting ^ actual 
lattice data for the effect of the hnite spatial box length L. 

^Here we assume that the finite volume correction is applied before the continuum and chiral extrapolations, 
thus interchanging steps (i) and (ii) of Sect. 1. In practice such a change is helpful, since otherwise the volumes 
would have to be matched to perform a well-defined continuum extrapolation and this would require a priori 
knowledge of the lattice spacing that will come out of the simulation. 
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Naive reasoning suggests that - because finite volume effects are due to the pion cloud 
around a particle and thus to pure IR physics - such shifts will be rather insensitive to the 
UV properties of the theory. This is what one expects to hold as long as the cut-off is large 
compared to the scale of chiral symmetry breaking, Axsb — 1 GeV. With a lattice regularization 
all momenta are cut off at vr/a, and with a standard lattice spacing a ~ 0.1 fm the resulting 
scale ~6 GeV is indeed much bigger than Axsb- 

This intuitive argument can be refined in two ways. The first option is to invoke an extension 
of GhPT designed to take care of the effects of the finite lattice spacing a. Of course, the details 
of this theory need to be tailored to the action used, but generically the new Lagrangian follows 
from the old one by replacing the low energy constants, e.g. £3 ^ + const W 3 - see Ref. jH] 

for a recent review. In consequence, the 0 {p'^) formula (jH)) takes the form 


M^{a,L) ^ M ^l-jX ^3 + const 103 - log(x) + 


( 88 ) 


where x = M'^/{4:nF)‘^, M = y/2Bm and £3 = log(A|/(47rF)^). In other words the very effect of 
such an extension concerns the particle mass in infinite volume, the fractional finite-size effect 
gets modified by a-effects only at O(p^) in the chiral counting, viz. 


Mjr(a, L) = Mt^^o) 


(1 + + Oix^)) . 


(89) 


The second option is to compare the generic one loop finite volume shift in the continuum 
to a version in which the pion propagator is discretized in the simplest ® possible way 


^i(M,L, 0 ) 


gi{M,L,a) 



1 
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1 

+ Po + 
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(90) 

(91) 


[the definition Si sin^(|pi) has been used and the finite sum runs over pi = ^Ui with 

Uj G {0, — and N = L/a an integer] and verify that the difference is small compared to 

the shift itself, i.e. \gi{M, L, 0)—gi{M, L, a)| <^gi{M, L, 0) for standard values of M, L, a. With 


gi{M, L, 0) 


gi{M,L,a) 



1 

lOvr^t^ 


^ m(n) 

n>l 


1 

dvr 2 


y~] m(n) 

n>l 


VnL 



r 1 

n=0 


At , , 7rn, \ 


- Jo(2t/a^) l3\ 

. 0 . j ^/A^^t 


(92) 

(93) 


[cf. Tab. [U for m{n)] and M = 300MeV, L = 2fm, a = 0.1 fm one finds {A 7 r/Mygi{M,L, 0 ) = 
0.4374 and {An/MfigfiM, L, a) =0.4400, thus a difference of less than a percent in a quantity 
designed to correct actual data by - at most - a few percent. It is hence sufficient to calculate 
the finite volume effects in continuum GhPT. This conclusion was also reached in Ref. [S]. 

®By considering mm we do not indicate that this would yield a better estimate of the finite size effects 
than the continuum form innina. The actual discretization of quarks and gluons will not lead to a simple pion 
propagator, but we are interested in the absolute difference \gi{M, L,0) — gi{M, L,a)\, since it is expected to 
correctly indicate the order of magnitude of discretization effects in the finite volume shift of actual data. We 
stress that our discretization is similar in spirit, but not identical to the one used in Ref. |42|. 
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C Effects due to kaon and eta loops 


Our formulae ((221121 take only the effects due to virtual pion loops into account. In other 
words, they neglect the contribution to Mp{L)—Mp and Fp{L)—Fp coming from kaon and eta 
loops “around the world”. In this appendix we want to discuss to which extent this is justihed. 

Consider a soon-to-be standard Nf =2 + 1 simulation with Mn- = 300MeV and the strange 
quark hxed at its physical value and Mk = 530MeV in consequence (see Fig.lH). With a box 
size L = 2 fm a hrst crude estimate says that the kaon loop effects will be down, relative to the 
pion loops, by a factor ^ 10 % correction on the fractional hnite volume 

effect is not exactly small. However, this correction should be compared to the absolute size of 
the effect and the statistical error and from Tabs. EHH and/or Figs. l21IKl it follows that it is safe 
to neglect such a correction at the permille level. Finally, we mention that one expects NNLO 
contributions for Rf„, to be of the same order of magnitude as for Rm.^ and this means that 
an additional pion loop could prove more important than a single kaon loop in hnite volume. 

We have verihed that kaon and eta loops “around the world” prove numerically insignihcant 
by comparing the pion-loop contributions to those of the kaons and etas in the full 1-loop 
expressions calculated in SU (3) ChPT in hnite volume. We hnd 


RMtt (-^tt) 


Rmk 


1 ^ , 1 . , , 

r 1 _ 1 _ 1 „ 1 

^Mr, = 2^Kgi{^K) - 

(A-tt) T g^x5'i(-^i^) T 


where Ap and gi have been dehned in ina and ca, respectively, and 


(94) 

(95) 

(96) 


Rf., = - ■^iKgii.^K) (97) 

3 _ 3 _ 3 _ 

^Fk g'C ttS'i (Att) g'C»?l?i (-^i?) (98) 

RFr, = -^^Kgii^x) ■ (99) 


These formulae deserve a few comments. First, a few elementary checks: Our (ED, (ED and 
agree with AMp'/Mp-, AFp'/Fp' as given in [Till IT7j and both the Rm and the 

Rf become degenerate in the SU{3) limit (m„ = rrid = mg), and furthermore, these degenerate 
expressions agree with the result by Gasser and Leutwyler |2j , eqns. (|HD and ©, specihed to 
Nf =3. Second, Rmk Rf^ have no ^(Aj^) (in other words only kaon- and eta-loops contribute 
at one-loop order) and in Rm^ the one-pion-loop contribution is suppressed by an extra factor 
Therefore, one expects Mk{L) — M^, Mrj{L) — and Fr^{L) — F^j to be small, 
and the numerical investigation in Sect. 6 specihes to which extent this is true. The numerical 
discussion also shows that, for a substantial range of pion masses, we have M^(L) < M^, 
in (apparent) contradiction to our statement in Sect. 2 that hnite volume effects will lift the 
masses of the pseudo Goldstone bosons. However, this just indicates that the general rule 
may be overwhelmed by SU{3) breaking effects. As our formulae show, the SU{3) breaking 
effects are accidentally dominating (at one-loop level) only in the Rmp and not in the Rfp, i-e. 
Rfp <0 for P = 7r,iF, 7]. Finally, we wish to elaborate on a point already raised in the work 
by Becirevic and Villadoro PI- The main message of the one-loop formulae (jHl-ESD is that 
hnite volume ehects and chiral logs are intimately related. In this particular case, where there 
is only a tadpole contribution, the hnite volume shift follows from the quark mass dependence 




by the simple substitution \og{Mp/jj?) \og{Mp/jj?) + gi{Xp). In practice, this means that 
one cannot extract “chiral logs” and the pertinent low energy constants without controlling the 
hnite volume effects. 
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